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Observation of Cold Collisions between Trapped Ions and Trapped Atoms 
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We study cold collisions between trapped ions and trapped atoms in the semiclassical (Langevin) 
regime. Using Yb + ions confined in a Paul trap and Yb atoms in a magneto-optical trap, we 
investigate charge- exchange collisions of several isotopes over three decades of collision energies 
down to 3 fieV (fcs x 35 mK). The minimum measured rate coefficient of 6 x 10 -10 cm 3 s" 1 is in 
good agreement with that derived from a Langevin model for an atomic polarizability of 143 a.u. 

PACS numbers: 34.70.+e, 31.15.ap, 32.10.Hq 



Studies of cold collisions between trapped neutral 
atoms have revealed a plethora of fascinating quan- 
tum phenomena, including Wigner threshold laws [l|, 
magnetically tunable Feshbach resonances [2] , controlled 
molecule formation 0], and the suppression of individ- 
ual scattering channels Collisions between trapped 
ions, on the other hand, are featureless, since the strong 
long-range repulsive Coulomb interaction prevents the 
ions from approaching each other. Collisions between 
ions and neutral atoms 0, @, 0, Q fall into an intermedi- 
ate regime where an attractive long-range r -4 potential 
leads to semiclassical behavior for a wide range of colli- 
sion energies, but where quantum phenomena dominate 
at very low energies. Cold ion-atom collisions have been 
proposed as a means to implement quantum gates |9[ , to 
cool atoms or molecules pjj, [12| lacking closed op- 

tical transitions, to bind small Bose-Einstein condensates 
to an ion [13], or to demonstrate novel charge- transport 
dynamics [141 ]. 

As a function of collision energy, ion-atom collisions ex- 
hibit three distinct regimes @, Q : a high-energy classical 
(hot) regime with a logarithmic dependence of cross sec- 
tion a on energy E, a wide semiclassical Langevin (cold) 
regime with a power-law dependence cr(E) oc E~ x l 2 
for charge-exchange and momentum transfer processes 
0], and a quantum (ultracold) regime where con- 
tributions from individual partial waves can be distin- 
guished. However, given the large forces on ions pro- 
duced by small stray electric fields, it has been difficult to 
reach experimentally even the semiclassical regime. Refs. 
[HI Gil , studying charge-exchange at E ~ 100 meV, re- 
port the only observations of Langevin- type ion-atom col- 
lisions. In ion- molecule systems, experimental signatures 
of Langevin collisions have been seen at high tempera- 
ture [13, III Ell, and recently also at 1 K (80 /xeV) 0. 
In all previous work, at most one of the collision partners 
was trapped. 

In this Letter, we study collisions between indepen- 
dently trapped, laser-cooled ions and atoms down to 
unprecedented low energy (3 /ieV) in the semiclassical 
collision regime. Using a double-trap system [2l|, we 



investigate resonant charge-exchange collisions for dif- 



ferent Yb + +Yb isotope combinations and find agree- 
ment with the Langevin model to within a factor of 
two over three decades of energy [H, @. The highest 
energy 4 meV = ks x 45 K corresponds to the transi- 
tion to the classical regime 22|, |23|, while at the low- 
est energy 3 /ieV = fc# x 35 mK, where approximately 
40 partial waves contribute to the cross section, isotope 
shifts should become relevant. The lower limit on colli- 
sion energy is set by our ability to detect and minimize 
the micromotion of a single ion in the Paul trap. 

The long-range interaction potential between a singly 
charged ion and a neutral atom is the energy of the in- 
duced atomic dipole in the ion's electric field, given by 
V(r) = — C4/(2r 4 ), where C4 = aq 2 / (Aireo) 2 is propor- 
tional to the atomic polarizability a, and q is the elec- 
tron charge. For given collision energy E in the center- 
of-mass frame, there exists a critical impact parameter 
b c = (2C4/E) 1 / 4 that separates two types of colli- 
sions: those with impact parameter b < b c that result 
in inward-spiraling orbits of radius approaching zero, and 
those with b > b c that never cross the angular- momentum 
barrier [H|, H3] . 

For a resonant charge-exchange process A + + A — > 
A + A + of an ion with its parent atom, a semiclassical 
cross section <j ce can be simply derived, provided that 
b c is large compared to the range of the molecular po- 
tential: for collisions with b > b c the electron should 
remain bound to the incoming atom, while for inward- 
spiraling collisions with b < b c the electron should, 
with equal probability, attach to either nucleus. It fol- 
lows that the resonant charge-exchange cross section is 
^"ce = pl/^, where <tl = 7rb 2 is the Langevin cross sec- 
tion. The corresponding semiclassical rate coefficient 
Kce — &ceV — Try/ C4J where v is the relative veloc- 
ity and \i the reduced mass, is independent of energy. 

At high energies the model becomes invalid when b c 
becomes so small that charge exchange outside the cen- 
trifugal barrier starts to play a rolejfor Yb, this occurs at 
energies E > 10 meV= k B x 120 K 0,111123. At very low 
energies the semiclassical Langevin model ceases to be 
valid when the s- wave scattering limit is reached, which 
happens near E s ~ h 4 /(2fi 2 C4 z ) = 4 peV = fc# x 50 nK. 
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Figure 1: Setup for trapping neutral (MOT) and singly 
charged (surface planar Paul trap) Yb in the same spatial 
volume. Stray dc field compensation and trapping along z for 
the Paul trap is provided by dc electrodes. Imaging CCDs 
point along (1, 0, 1) and (0, 1, 1) axes, (color online) 



Quantum simulations [6[ of this system [22[ confirm 
the semiclassical Langevin model in the indicated en- 
ergy range. For collisions between different isotopes, the 
Langevin expression should be modified at low collision 
energies when they become comparable to the small dif- 
ference in binding energy of the electron to the two nuclei, 
i.e. the isotope shift of the ionization potential. In this 
case we expect (exo-) endo-energetic charge exchange col- 
lisions to be (enhanced) suppressed [25|, [26[ . 

We use a magneto-optical trap (MOT) in combination 
with a Paul trap, as originally proposed by W. Smith [ll| . 
The setup, shown in Fig. [TJ is described in detail else- 



where [21[. In the present work, 172 Yb, 174 Yb, or 171 Yb 



atoms are selectively loaded from an atomic beam into a 
MOT by tuning the laser frequency near the 1 So -^ x Pi 
transition at A a = 398.8 nm. Typically, the MOT is 
operated at a 75 G/cm magnetic- field gradient and con- 
tains 3xl0 5 atoms at a peak density of 2 x 10 8 cm -3 and 
temperature of 700 fiK as determined by a time-of-flight 
measurement. 

Cold ions are produced by nonresonant photoioniza- 
tion from the excited x Pi state of the MOT with 370- 
nm light from a secmiconductor laser |2l|. The ion trap 
is a surface-electrode Paul trap printed on a vacuum- 
compatible substrate 21], and is typically operated at 
1.4 MHz to create a 0.3 eV deep pseudopotential trap 
3.6 mm above the trap surface with a secular frequency 
of 67 kHz. Ion trap populations can be adjusted between 
a single and 10 4 ions by varying the trap loading time. 

Two beams from the same 370 nm laser used for pho- 
toionization provide Doppler cooling of the ions on the 
2 Si/2 ^ 2 Pi/2 transition along all three principal trap 
axes. Ions that decay to a met ast able D state are re- 
pumped [27] using a laser operating at 935 nm [2l|. We 
detect the ion population by monitoring trap fluorescence 
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Figure 2: Typical 172 Yb + ion- crystal fluorescence decay and 
exponential fits for no (red circles), moderate (blue squares), 
and high (green triangles) 174 Yb atomic density at the ion 
trap site. Inset: 172 Yb + (dashed, gold) and 174 Yb + (solid, 
purple) population evolution in the presence of 174 Yb. Both 
traces are normalized by the same peak value, (color online) 



at 370 nm with a photomultiplier tube. A single cold 
trapped ion produces 5 kcounts/s. 

Collisions that only change the particles' energy and 
momentum are difficult to observe in our setup due to 
the continuous laser cooling. On the other hand, charge- 
exchange collisions between different isotopes a Yb + + 
^Yb -^ a Yb + ^Yb + are easy to observe using isotope- 
selective ion fluorescence: we first load the ion trap from 
the MOT with isotope a Yb+, change the MOT isotope 
to ^Yb by adjusting the frequency of the 399-nm laser, 
and then monitor the decay of the a Yb + ion popula- 
tion through the decay of the 370-nm fluorescence. In 
order to prevent photoionization and direct loading of 
the ion trap with ^Yb + , we modulate the MOT and ion 
light out of phase to ensure that the MOT contains no 
excited atoms when the 370-nm light is present. With- 
out the MOT, the ion trap loss is exponential with a 
typical lifetime To = 400 s (Fig. [2j red circles) for ion 
crystals, presumably due to collisions with background 
gas atoms. In the presence of the MOT, r is substan- 
tially shortened, with a shorter lifetime for higher MOT 
density. The simple exponential decay over a decade in 
ion number indicates that the loss involves a single ion, 
rather than collisions between ions. We also measure the 
trap fluorescence for identical ion and MOT isotopes. In 
the latter case, we observe a small initial decay, presum- 
ably due to a small amount of heating of the outermost 
ions, but no exponential loss (Fig. [2] inset). This confirms 
that the measured trap loss for different ion and MOT 
isotopes is indeed due to charge-exchange collisions, and 
not to heating in ion-atom collisions that interrupt the 
ion's micromotion [28| . 

To accurately determine the charge-exchange rate coef- 
ficient, we vary the atomic density at the ion trap location 
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Figure 3: (a) (1,0, 1) camera image of the ion crystal (blue) 
and cross-section showing highly non-Gaussian shape of crys- 
tal (red), (b) (0,1,1) camera image of the ion crystal, (c) 
Typical low-overlap setting between MOT (colored contours) 
and 1/e 2 contour of ions (white), (d) Same as (c) but for a 
higher overlap setting, (color online) 



by moving the zero of the MOT magnetic quadrupole 
field with a bias field. We determine the local atomic 
density at the ions' location by taking images of both 
the MOT cloud and the ions on charge-coupled de- 
vice (CCD) cameras along two different directions (see 
Fig. [3]), with the total atom number calibrated by means 
of a resonant-absorption measurement. We then calcu- 
late the average atomic density experienced by the ions 
(n) = J p(r)n(r)dr, where p(r) is normalized ion distri- 
bution and n (r) is the local atomic density. The inelastic 
rate coefficient K is obtained as the slope of a linear fit 
of the observed ion decay rate constant Y = 1/r vs. (n). 

We note that in most situations the ions' kinetic energy 
is determined by their micromotion throughout the rf 
cycle [29]. Any ion displacement from the zero of the 
oscillating electric field, be it due to ion crystal size or 
a dc electric field, results in micromotion whose energy 
easily exceeds that of the thermal motion in the secular 
potential of the Paul trap. To investigate the dependence 
of the rate coefficient K on collision energy E, we vary 
the latter by loading different numbers of ions into the 
trap or by intentionally offsetting the ion crystal from the 
zero of the rf quadrupole field with a dc electric field. 

The rate coefficient K(E) is shown in Fig. [U In the 
region E > 30 fieV the ions' energy is calculated from 
the observed Doppler broadening. For well-compensated 
traps containing small crystals or single ions, where the 
Doppler broadening is smaller than the natural linewidth, 
we determine the average micromotion energy from the 
observed correlation between rf drive signal and ion fluo- 
rescence [29|. In the latter case, the signal-to- noise ratio 
is small, resulting in larger energy uncertainty. For all en- 
ergies investigated here, the atoms' contribution to the 
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Figure 4: Charge- exchange rate coefficient K as a func- 
tion of collision energy in the center-of-mass frame. Red 
circles, green, and blue diamonds represent 172 Yb + + 174 Yb, 
172 Yb + + m Yb, and 174 Yb + + 172 Yb, respectively. The solid 
line indicates the theoretical Langevin rate coefficient K ce 
[22l ], the dashed line assumes a contribution from the the P- 
state polarizability [3(| (see text). The black arrow indicates 
the ionization isotope shift between 174 Yb and 172 Yb. (color 
online) 



collision energy is negligible. 

The data point at the lowest energy E = 3.1 /ieV, be- 
ing an average of 22 single-ion measurements, has large 
uncertainty in K. It also has large energy uncertainty 
because the micromotion is sufficiently well compensated 
for the fluorescence correlation signal to be dominated by 
noise. The other data points show observable Doppler 
broadening and were measured with up to a few thou- 
sand ions, resulting in smaller statistical uncertainties. 
We estimate systematic uncertainties to be a factor of 
two in K due to the difficulty of absolute MOT density 
calibration, and ±50% in E, due to the non-thermal en- 
ergy distribution of the micromotion. 

The ab initio calculated value for the polarizability 
a s =143 a.u. = e -2.66 x 1(T 28 m 3 [22] for the atomic % 
ground state yields K ce = 5.8xl0 -10 cm 3 s _1 for ground- 
state collisions (solid line in Fig. [4]). Our experimentally 
measured value of K = 6 x 10 -10 cm 3 s _1 in the middle of 
the semiclassical region [22| around E ~100 fieV is thus 
in good agreement with the Langevin model. The data 
point at the highest collision energy E = 4 meV coincides 
with the transition to the classical region, where the rate 
coefficient increases with collision energv 0,(231. Doubly 
excited collisions Yb + *+Yb* cannot occur as the excita- 
tion light for ions and atoms is modulated out of phase. 
For Yb + *+Yb collisions, the Langevin cross section, de- 
pending only on the ion's charge and the atom's po- 
larizability, is unchanged, and the charge-exchange pro- 
cess remains resonant, yielding the same value of <j ce . 
Yb + +Yb* collisions cannot contribute for our collision 



energies since the atomic transition is tuned out of res- 
onance with the MOT light at an ion-atom distance 
Ro ~ 30 nm larger than b c = 2 — 12 nm, and modeling 
of collision trajectories shows that the time required to 
move between r = Rq and r = b c exceeds several excited- 
state lifetimes. For illustration purposes, we have indi- 
cated in Fig.[4]the quite small change in K ce if Yb + +Yb* 
collisions with o;p = 500 a.u. for the 1 Pi state [30] were 
to contribute at our measured MOT excited-state frac- 
tion. 

The theoretical value K ce = ctl^/2, corresponding to 
equal binding probability of the electron to the two nu- 
clei, should apply when the collision energy far exceeds 
the isotope shift of the ionization potential. The lat- 
ter can be estimated from spectroscopic data of a tran- 
sition to a state with low electron probability density 
at the nucleus, such as 4/ 14 6s 2 -> 4/ 14 6sl0<i It 
follows that the 172 Yb++ 174 Yb 172 Yb+ 174 Yb+ reac- 
tion should be exothermic and release AE = 2.9 fieV = 
h x 0.7 GHz, and other isotope combinations give sim- 
ilar values. At the substantially larger collision energy 
E = 0.21 meV ^> AE we have investigated various iso- 
tope combinations, 172 Yb++ 174 Yb, 174 Yb++ 172 Yb, and 
174 Yb + + 171 Yb, and we indeed find that they they all dis- 
play the same rate coefficient (see Fig. |4]). For the data 
point 172 Yb++ 174 Yb 172 Yb+ 174 Yb+ at the lowest en- 
ergy E = 3.1 /ieV ~ AE, we speculate, based on calcula- 
tions for similar processes in H [H, [26[ , that the exother- 
mic collision rate should exceed the resonant Langevin 
value, which is (weakly) supported by our data. 

While it may be impossible to compensate stray fields 
well enough to reach the s-wave scattering limit E s = 
4 peV, Feshbach [8] and other collision resonances may 
be observable well above E s 0, [22[ . All the Yb isotopes 
used here have been cooled to quantum degeneracy in 
an optical dipole trap [H, [33[, which would allow the 
investigation of collision processes between an ion and 
a Bose-Einstein condensate or Fermi gas. Alternatively, 
to avoid the large resonant charge-exchange cross section 
observed here, a different species such as Rb could be 
used for sympathetic cooling of ions @, EqI , or for study- 
ing ion impurities in a Bose-Einstein condensate [T3 |. 
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